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The envelope glycoproteins S1 and S2 of severe acute respiratory syndrome coronavirus (SARS-CoV) 
mediate viral entry by conformational change from a prefusion state to a postfusion state that enables 
fusion of the viral and target membranes. In this work we present the characterization of the dynamic 
properties of the SARS-CoV S2-HR2 domain (residues 1141-1193 of S) in the prefusion and newly discov- 
ered transition states by NMR !°N relaxation studies. The dynamic properties of the different states, 
which are stabilized under different experimental conditions, extend the current model of viral mem- 
brane fusion and give insight into the design of structure-based antagonists of SARS-CoV in particular, 
as well as other enveloped viruses such as HIV. 


© 2009 Elsevier Inc. All rights reserved. 


SARS-CoV! (severe acute respiratory syndrome coronavirus) en- 
ters target cells by first binding to the cell membrane and subse- 
quently inducing fusion of the viral and cellular membranes [1,2]. 
The envelope proteins S1 and S2 mediate the binding and fusion 
steps, respectively [3], in a manner that is thought to be similar to 
the well-characterized envelope proteins of HIV and influenza 
[4,5]. During the fusion step, SARS-CoV S2, as well as its homologs 
in HIV and influenza, is thought to undergo a series of conforma- 
tional changes from a prefusion state to a postfusion state [1-5]. 
SARS-CoV S2 possesses 2 heptad repeat domains termed HR1 and 
HR2, which have a high propensity to form helix. The high-resolu- 
tion structure of the SARS-CoV S2 in the postfusion state has been 
determined [6-8]. In the postfusion state of SARS-CoV S2, as well 
as the postfusion form of the analogous proteins in Ebola, HIV, influ- 
enza and parainfluenza [9-12], the HR1 domains self associate as 
parallel trimers and the HR2 domains associate with the HR1 in an 
anti-parallel orientation to form the “6 helix bundle”. More recently 
structures of the prefusion states of SARS-CoV S2-HR1 [8] and S2- 
HR2 [13] have been determined. Other examples of prefusion struc- 
tures include influenza HA2 [14] and parainfluenza virus 5 F [15]. 
Support for the presence of the prefusion conformation in vivo arises 
from the anti-viral activities of peptides, corresponding to the HR1 
and HR2 sequences, that apparently inhibit virus entry by binding 
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to the prefusion conformation and thereby prohibit the transforma- 
tion to the postfusion form |16-18]. A peptide termed T20 or fuzeon, 
which stems from the HIV gp41 HR2 sequence, is currently used in 
humans [19]. 

A central feature of the viral entry model is the transformation 
from the prefusion state to the postfusion state. Recently, S2-HR2 
has been shown to be in a dynamic equilibrium between a struc- 
tured coiled-coil trimer, which is thought to represent the prefu- 
sion state, and unstructured monomer, which is thought to 
represent a metastable transition state [20]. In the present work, 
we present the detailed characterization of the dynamic properties 
of SARS-CoV S2-HR2 in the prefusion and transition states using 
NMR relaxation. We will discuss the importance of considering 
the S2-HR2 dynamic properties in the viral entry model, as well 
as in the design of future envelope-based anti-virals. 

Recombinant SARS-CoV S2-HR2, consisting of residues 1141- 
1193 of S, was prepared as previously described [13]. NMR samples 
were prepared by ultrafiltration (YM3, Amicon, Billerica, MA) using 
the appropriate buffer (10 mM PO,/pH 7.0 or 10 mM PO,/pH 7.0, 
30% TFE-d3, Cambridge Isotope Laboratories, Andover, MA). The fi- 
nal concentration of S2-HR2 was 400 uM (Ry, of the prefusion state 
was also determined at 1 mM protein concentration). NMR spectra 
were recorded on a Bruker AVANCE 800 MHz spectrometer 
equipped with a cryogenic triple resonance probe. Spectra were 
processed by NmrPipe and visualized with NmrDraw [21]. The 
backbone assignments at 25°C in the presence and absence of 
30% TFE have been reported elsewhere [13,20]. The backbone 
assignments of S2-HR2 at 45 °C were determined from a series of 
the standard triple resonance experiments [13]. The '°N R,(1/T,), 
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Rxy(1/T2) and HNOE values were measured by standard pulse se- 
quences [22]. For determination of R,, the delay values were set 
to: 30, 60, 120, 240, 72 and 1250 ms (prefusion state); 50, 100, 
200, 400, 800 and 1200 ms (transition state). For determination 
Of Rxy, the delay values were set to: 8, 16, 32, 64, 96, 128, 256 
and 512 ms (prefusion state); 10, 20, 40, 80, 128 and 256 ms (tran- 
sition state). The decays of cross-peak intensities with time in the 
'N R, and Ry experiments were fit to a single exponential by a 
non-linear least-squares fit with fitting errors typically less than 
5% (Kaleidagraph 3.08, Synergy Software). No significant difference 
for the prefusion R,, was observed for two different concentrations, 
400 1M and 1 mM, suggesting that R,, was not dependent on con- 
centration. HNOE values for each given residue were calculated as 
the intensity ratio (I/Ip) of the '°N-'H correlation peak in the pres- 
ence (J) and absence (Jj) of proton saturation during the relaxation 
delay of 5s. Errors in the HNOE were estimated to be +0.10 from 
the S/N ratios. For reduced spectral density mapping, the spectral 
density functions were described as: 


J(O) = (6Ry — 3R, — 2.720nn)/(3d° + 4c?) 
J(@n) = (4R; — 5Ony)/(3d? + 4c?) 


J(0.87@y) = 40ny/(5d7) 


where ony = RANOE — 1)yn/Vuj d = (Moh ynyu)/(82217); c= @NAo/31! 
2.» and yy are the gyromagnetic ratios of '°N and 'H, respectively; 
[lo is the free space permeability; h is Planck’s constant divided by 
27; ris the length of the 'H-!°N bond, 1.04 A; my and wy are the 
Larmor frequencies of '°N and 'H, respectively; Ao is the '°N CSA 
of —162 ppm [23]. Errors of the spectral density functions were 
determined by propagating the uncertainties of the individual 
relaxation parameters. 

S2-HR2 exists in different conformational states, which are 
dependent upon the experimental conditions. For example, at 
25 °C in the presence of 30% TFE, S2-HR2 is a stable coiled-coil tri- 
mer, which represents the prefusion state [13]. On the other hand, 
in the absence of co-solvent, at 25 °C S2-HR2 is in equilibrium be- 
tween unstructured monomer and a structured trimer, which is 
thought to represent an equilibrium between the prefusion and 
transition states [20]. Moreover, in the absence of co-solvent, S2- 
HR2 exhibits a melting temperature of ~35 °C and is completely 
unstructured at 45 °C [20]. To gain insight into the dynamic prop- 
erties of the prefusion and transition states, we characterized the 
dynamic properties of S2-HR2 under two experimental conditions: 
10mM NaPO,/pH 7.0, 30% TFE at 25°C (prefusion state) and 
10 mM NaPO,/pH 7.0 at 45 °C (transition state). NMR relaxation 
is correlated to the global correlation time t,, which can be esti- 
mated from the Stokes equation (t, = (471r°7)/(3k,T), where r= ra- 
dius, 7 = viscosity, k, = Boltzman’s constant and T= temperature). 
Consequently, the interpretation of the NMR relaxation properties 
of S2-HR2 must take into consideration the temperature, solvent 
viscosity, and protein molecular volume, of S2-HR2 as summarized 
in Table 1. For example, the estimated Tt, of the transition state (i.e. 
unstructured monomer) is predicted to be a factor of ~4 less than 
that of the prefusion state (i.e. structured trimer in the presence of 
co-solvent). 

The dynamic properties of S2-HR2 under two different experi- 
mental conditions were characterized by determination of hetero- 
nuclear relaxation rates of the NH, which have been widely used to 
investigate the backbone dynamics of proteins and the role of 
dynamics in protein function [26,27]. The '°N R, and R,, relaxation 
rates and the 'H-!°N NOE (HNOE) at 800 MHz of S2-HR2 are 
shown in Fig. 1 (cf. Supplementary Tables 1 and 2). The average 
R,, Ryy and the HNOE values of the prefusion S2-HR2 helix (residues 
17-47 of the construct, residues 1155-1185 of S) are summarized 


Table 1 
Predicted molecular dynamic properties of SARS-CoV S2-HR2 (residues 1-55 of the 
construct, residues of 1141-1193 of S). 


State? T (°C) Mw (kDa) n (cP) CARER 
Transition 45 6.3° 0.61 0.23 
Prefusion 25 18.6¢ 1.63 1.00 


* Transition = 10 mM NaPO,/pH 7.0 at 45 °C; prefusion = 10 mM NaPO,/pH 7.0, 
30% TFE at 25 °C. 

> (te =(Anr?n)/(3kpT), where r = radius, 7 = viscosity, kp, = Boltzman’s constant and 
T = temperature) and Tye = T. Of the structured trimer (i.e. at 10 mM NaPO,/pH 7.0, 
30% TFE at 25 °C). Rmonomer is taken to be 23.1 A, which is based on the experimental 
value for protein G, a protein with an identical number of residues, under dena- 
turing conditions [24]. Rtrimer is taken to be 26.7 A, as calculated by CNS [25] using 
the coordinates for the prefusion form [13]. 

© Based on monomer-trimer equilibrium in fast exchange at 25 °C (70% mono- 
mer/30% trimer) and the lack of structure at 45 °C (100% monomer) [20]. 

¢ Based on the solution structure and analytical centrifugation studies [13]. 


in Table 2. The values for the prefusion state suggest a helix exhib- 
iting a relatively rigid region with increased flexibility at its C-ter- 
minal end, and flexible N- and C-terminal regions. The relaxation 
values for the transition state are consistent with the faster molec- 
ular tumbling of a smaller molecule under conditions of higher 
temperature and lower solvent viscosity (Table 1). The relative uni- 
formity and magnitude of the transition state relaxation parame- 
ters for residues 17-47 support the notion that under these 
conditions S2-HR2 is monomeric with little residual structure 
present. Nonetheless, the termini of the transition state are clearly 
more flexible on the ps—ns timescale as indicated by the negative 
HNOE. 

The reduced spectral density functions J(0), J(@y) and J(0.87@ ) 
are correlated with the magnitude of NH bond fluctuations at three 
different frequencies. In contrast to the widely used model free for- 
malism, spectral density functions make no assumptions for a mo- 
tional model with a uniform Tt, and thus spectral density functions 
are clearly most appropriate for the analysis of the dynamic prop- 
erties of partially unfolded and unstructured regions [23], which is 
clearly relevant in the case of S2-HR2. In Fig. 2, the J(0), J(@y) and 
J(0.87@) values are shown for the residues of S2-HR2 under the 
two experimental conditions. The average spectral density values 
of the prefusion S2-HR2 helix (residues 17-47 of the construct, res- 
idues 1155-1185 of S) are summarized in Table 2. In the prefusion 
state, higher J(0) values are found in the helix and lower J(O) are 
found at the termini, which is consistent with the presence of a 
coiled-coil core on the ns timescale. In the case of the transition 
state, lower and more uniform J(0) are exhibited, consistent with 
an unstructured monomer. The uniform dynamic profile of the 
S2-HR2 transition state resembles that of the urea-denatured pro- 
tein G [28]. Note that the ratio of the J(0) values can be used to esti- 
mate the ratio of the t. values. For example, the ratio of the helix 
region of the transition state /(0) to that of the prefusion J(O) is 
0.24, which is in very good agreement with that predicted in Table 
1 (the prefusion state tT, is predicted to be ~4~ that of the transi- 
tion state T,), especially considering the difficulty in estimating the 
molecular volume of the unstructured monomer. With respect to 
the high frequency motions reflected by J(@y) and J(0.87@wn), the 
prefusion state exhibits significantly lower ps motions in the helix 
with respect to the transition state. However, the C-terminal region 
of the prefusion helix, consisting of ~8 residues, exhibit increased 
dynamic motions similar to that observed for J(0). Interestingly, 
this region is not helical in the absence of co-solvent [20] and thus 
the increased dynamics may be reporting a region of inherent 
instability. In contrast, NMR dynamic studies have generally shown 
that helices exhibit uniform fast time scale dynamics with small 
increases in ns—ps motions at the terminal three residues of the he- 
lix [29]. A notable exception has been observed for the basic leu- 
cine zipper of GCN4, which is another coiled-coil helix [30]. 
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Fig. 1. '°N relaxation parameters of SARS-CoV S2-HR2 under two different 
experimental conditions: (a) '°N R,, (b) '°N Ry and (c) HNOE. Experimental 
conditions were 400 uM S2-HR2 in 10mM PO,/pH 7.0, 30% TFE at 25°C (filled 
circles) and 400 uM S2-HR2 in 10 mM PO,/pH 7.0 at 45 °C (open circles). 


Table 2 
Average relaxation parameters in the SARS-CoV S2-HR2 helix (residues 17-47 of the 


construct, residues 1155-1185 of S). 


Relaxation parameter Prefusion* Transition” 
Raises) 0.95 + 0.20 1.64 + 0.10 
Ry (s~') 23.5+4.7 6.33 + 1.39 
HNOE 0.66 + 0.07 0.08 + 0.16 
J(O) (ns/rad) 5.44 + 1.12 1.32 + 0.37 
J(@n) (ns/rad) 0.14 + 0.03 0.23 + 0.02 
J(0.87@y) (ns/rad) 0.0052 + 0.0019 0.023 + 0.003 


* Experimental conditions were 400 uM S2-HR2 in 10 mM NaPO,/pH 7.0, 30% TFE 


at 25°C. 


> Experimental conditions were 400 1M S2-HR2 in 10 mM NaPO,/pH 7.0 at 45 °C. 
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Fig. 2. Reduced spectral density mapping of SARS-CoV S2-HR2 under two different 
experimental conditions: (a) J(0), (b) J(@y) and (c) J(0.87@y). Experimental 
conditions were 400 uM S2-HR2 in 10mM PO,/pH 7.0, 30% TFE at 25 °C (filled 
circles) and 400 uM S2-HR2 in 10 mM PO,/pH 7.0 at 45 °C (open circles). 


In summary, the relaxation properties of S2-HR2 are consistent 
with the model of an equilibrium between a prefusion state with a 
coiled-coil trimer conformation and a transition state comprised of 
unstructured monomer [20]. In Fig. 3, such a model is presented 
with the structural and dynamic properties of the two states sum- 
marized. An important feature of this model is that the transition 
state enables the conformational changes of envelope necessary 
for SARS-CoV entry (i.e. interaction between HR2 and HR1 to form 
the 6 helix bundle). Based on the present characterization of the 
S2-HR2 dynamic properties, the transition state may accomplish 
this by efficient sampling of a large area of conformational space. 
Finally, we note that, due to the very different structural and dy- 
namic properties of the prefusion and transition states, consider- 
ation of the their structural and dynamic properties may be of 
importance in the design of HR2-based anti-virals for SARS-CoV 
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Fig. 3. Model of SARS-CoV S2-HR2 prefusion and transition states. The average relaxation parameters shown are those for the prefusion helix (residues 17-47 of the 
construct, residues 1155-1185 of S) and the analogous region of the transition state. Note that the prefusion and transition state values were acquired at different 
temperatures and solvent viscosities (cf. Table 1). The relative locations of the HR1 and transmembrane domains, which occur at the N- and C-termini of S2-HR2, respectively, 
are depicted by dotted lines. The structure of the prefusion state is taken from Hakansson-McReynolds et al. [13] with each subunit colored differently; the structure of the 
transition state was modeled by torsion angle simulated annealing of the S2-HR2 monomer in the absence of NMR restraints. The equilibrium constant shown was estimated 
from the observation that the '°C secondary chemical shifts of the S2-HR2 helix at 25 °C are ~30% of the magnitude observed in the presence of TFE [20]. 


in particular and perhaps analogous HR2-based anti-virals for 
other viruses such as HIV. 
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Appendix A. Supplementary data 


Two tables containing the relaxation data (R,, Ry, and HNOE) 
and the reduced spectral density values (J(0), J(@n) and 
J(0.87c};)) for SARS-CoV S2-HR2 under different experimental con- 
ditions. Supplementary data associated with this article can be 
found, in the online version, at doi:10.1016/j.jmr.2009.09.012. 
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